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Abstract 

A globally exponentially stable speed observer for mechanical systems was recently reported in the literature, under the 
assumptions of known (or no) Coulomb friction and no disturbances. In this note we propose and adaptive version of this 
observer, which is robust vis-a-vis constant disturbances. Moreover, we propose a new globally convergent speed observer 
that, besides rejecting the disturbances, estimates some unknown friction coefficients for a class of mechanical systems that 
contains several practical examples. 
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1 Introduction 

The design of speed observers for mechanical systems is a 
problem of great practical importance that has attracted 
the attention of researchers for over 25 years—the reader 
is referred to the recent books [1,3] for an exhaustive list 
of references. The first globally exponentially convergent 
speed observer for general, simple mechanical systems 
was recently reported in [2], where the Immersion and 
Invariance (I&I) techniques developed in [1] were used. 
Although the observer of [2] considers the case of sys¬ 
tems with non-holonomic constraints, it relies on the as¬ 
sumptions of no friction and the absence of disturbances. 
In [13] this speed observer was redesigned to accommo¬ 
date the presence of known Coulomb friction and it was 
used to design a uniformly globally exponentially stable 
tracking controller using only position feedback for me¬ 
chanical systems (without non-holonomic constraints). 

In this paper we propose two new robust velocity ob¬ 
servers for mechanical systems (without non holonomic 
constraints). First, we add an adaptation stage to the 
observer of [2] to reject constant input disturbances. Sec¬ 
ond, for mechanical systems with zero Riemann symbols 
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(ZRS) [11,14,16], we propose a new adaptive speed ob¬ 
server that, besides rejecting the disturbances, estimates 
some of the unknown friction coefficients. It should be 
noted that in this case there are products of unmea¬ 
surable states and unknown parameters, a situation for 
which very few results are available in the observer de¬ 
sign literature—even for the case of linear systems. A 
similar transformation has been presented in [4,5], where 
a coordinate transformation that removes the quadratic 
terms in velocity is found to solve challenging position- 
feedback tracking problems for surface ships and mobile 
robots. 

The paper is organized as follows. The two adaptive ob¬ 
servation problems addressed in the paper are presented 
in Section 2. The standing assumptions and a prelim¬ 
inary lemma are given in Section 3. In Section 4 we 
present an adaptive observer for systems with unknown 
friction and disturbances. In Section 5 the I&I observer 
of [13] is redesigned to accommodate the possible pres¬ 
ence of disturbances and known friction. Some physical 
examples are given in Section 6. The paper is wrapped- 
up with some future research in Section 7. 

Notation. To avoid cluttering the notation, throughout 
the paper n and a are generic positive constants. I n is the 
nxn identity matrix and O raxs is an n x s matrix of zeros, 
O n is an n-dimensional column vector of zeros. Given 
ai £ R, i € n := {1,..., n}, we denote with col(<+) the 
n-dimensional column vector with elements ai. For any 
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matrix A £ R nxn , {A)i £ R” denotes the *-th column, 
(. A) 1 the i-th row and ( A)ij the y'-tli element. That is, 
with e,i £ R”, i £ n, the Euclidean basis vectors, ( A)i := 
Aa, ( A) 1 ■= ejA and (A)ij := ejAej. For x £ R”, 
S £ R nxra , S = S T > 0, we denote the Euclidean norm 
\x\ 2 := x T x, and the weighted-norm ||a;||| := x T Sx. 
Given a function / : R™ —» R we define the differential 
operators 


V/:= 



V,,/ 



where Xi £ R p is an element of the vector x. For a map¬ 
ping g : R™ —> R m , its Jacobian matrix is defined as 


V.9 := 


(V.gi) T 

(V.g m ) T 


where gi : R ra —> R is the i-th element of g. 


2 FORMULATION OF TWO ROBUST SPEED 
OBSERVATION PROBLEMS 

In the paper we consider ?r-degrees of freedom, per¬ 
turbed, simple, mechanical systems described in port- 
Hamiltonian (pH) form [15] by 



VR(g, P ) + 


0 



0 

d 


(1) 


with total energy function H : R" x R" —>■ R 

H(q,p) = ^p r M-\q)p + V(q), (2) 

where q, p £ R” are the generalized positions and mo¬ 
menta, respectively, u £ R m is the control input, G : 
R” —> R raxm is the input matrix, the inertia matrix 
M : R n —>• R nxn verifies M(q) = M T (q) > 0 and 
V : R™ —> R is the potential energy function. As cus¬ 
tomary in the observer literature, it is assumed that the 
control signal u(t) is such that trajectories exist for all 
t > 0. 


The system is subject to two different perturbations. 

- Unknown constant disturbances d = col(di) £ R”. 

- Coulomb friction captured by 


IK = diag{n,r 2 , .,r„} £ R nxn , (3) 


with unknown r* > 0, i £ n. 


The problem is to design a globally convergent robust 
adaptive observer for the momenta p. The main contri¬ 
butions of the paper are the following. 

(i) For systems with ZRS design a new observer that is 
globally convergent in spite of the presence of the dis¬ 
turbance d and some unknown friction coefficients . 

(ii) If the friction is known robustify the observer of [2] 
to reject the disturbance d. 


Remark 2.1 The qualifier “some" in item (i) is essen¬ 
tial because, as will become clear later, except for the 
case of constant inertia matrix, we will not be able to 
consider the presence of unknown friction in all general¬ 
ized coordinates. 

Remark 2.2 Notice that the objective is to observe 
only the momenta (equivalently, the velocity) not to 
ensure consistent estimation of the parameters d and 
r := col(rj). As is well known in the identification liter¬ 
ature, a necessary condition for parameter convergence 
is that the signals satisfy a persistency of excitation con¬ 
dition [10]. In this respect, we notice that the system (1), 
(2), which can be represented in the state space form 


q = v 

v = U(q,v) — 9\v + G(q)u + d 


r = 0 
d = 0, 


for some U : R" x R n —> R™, with measurement q does 
not satisfy the observability rank condition [7] at zero 
velocity, hampering the observation of the parameters 

rB 

Remark 2.3 See Remark 1 in [12] for a physical inter¬ 
pretation of the disturbances d that, we underscore, en¬ 
ter at the level of the momenta. 


3 A Suitable pH Representation 

As shown in [16], the change of coordinates 

(<l,T T (q) p), 

where T : R™ —> R nxra is a full rank factorization of the 
inverse inertia matrix, that is, 

M-\q) = T(q)T T (q), (4) 


1 The authors thank Prof. Witold Respondek for this in¬ 

sightful remark. 
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transforms (1) into 


(iii) There exists a mapping Q : R" —> R" such that 


q 


i 

o 

1_ 

,p. 


J(q,p) - R(<i) _ 



0 


0 

+ 

T T (q)G(q)_ 

u + 

T T {q)d_ 


with new Hamiltonian W : R ra x R ra —>• R 
W{q,p) = i|p| 2 + V(q), 

the jk element of the skew-symmetric matrix J : R™ x 
R” —> R nx " given by 

J jk (q,p) = -p T [(T) j ,(T) k ], (6) 

with [•,•] the standard Lie bracket [14] and the trans¬ 
formed friction matrix 

R(q):=T T (q)XT(q)> 0. (7) 

Remark 3.1 One possible choice of the factorization 
(4) is the Cholesky factorization [6,8]. But, as will be¬ 
come clear below, other choices may prove more suitable 
for the solution of the problem. 

4 ROBUST OBSERVER FOR A CLASS SYS¬ 
TEMS WITH FRICTION AND DISTUR¬ 
BANCES 

In this section we solve the problem of robust observation 
of momenta in the presence of disturbances and some 
unknown friction coefficients for a class of mechanical 
systems. 


VQ(g) = T~\q). (9) 

Remark 4.1 Mechanical systems verifying Assump¬ 
tion 4.1 have been extensively studied in analytical 
mechanics and have a deep geometric significance— 
stemming from Theorem 2.36 in [11]. They belong to 
the class of systems that are partially linearizing via 
change of coordinates studied in [16]—see that paper 
for some additional references. 

4-2 Assumptions on friction 

To design our robust adaptive observer, besides Assump¬ 
tion 4.1, a restriction on the friction coefficients is im¬ 
posed. Namely, we assume that there are s, with s < n, 
unknown coefficient and decompose the friction matrix 
91 (3) as 

91 = 91 k + 91 u 

where 91k, 91 u are nxn diagonal matrices containing the 
known and the unknown friction coefficients respectively. 
As a working example consider the case n = 3 and s = 2 
with 

91 k = diag{0, r 2 , 0}, 9f u = diag{n, 0, r 3 } 

Similarly, with an obvious definition, we decompose the 
transformed friction matrix (7) into 

R(q) = Rk(q) + Ru{q)- 

To streamline the presentation all friction coefficients are 
grouped in a vector r = col(rj) £ R™ with the unknown 
and known coefficients in vectors r u £ R s andrfc £ R n-S , 
respectively. Thus, for our working example we have 

r = col(ri,r 2 ,r 3 ), r k = r 2 , r u = col(ri,r 3 ). 


4-1 Assumption on M(q) and a preliminary lemma 


Assumption 4.1 M l (q) admits a factorization (4) 
with a factor T(q) verifying 




0, i,j £ n. 


( 8 ) 


We also define a set of integers k C n that contains the 
indices of the unknown coefficients of r, which in the 
example is n = {1, 3}. 

Finally, we define a matrix C £ R nxs such that 

C T r = r u . (10) 


Instrumental for the developments of this paper is the 
following result, whose proof may be found in [16]. 

Lemma 4.1 The following statements are equivalent: 


Clearly, the matrix C verifies: 

• rank{C} = s. 

• For j £ k, ( C)j = e Kj . 


(i) M(q) satisfies Assumption 4.1. 

(ii) The Riemann symbols of M(q) are all zero PI 


2 See equations (6) and (7) of [16] for the definition of these 
symbols. 


In our example 


C = 


1 0 
0 0 
0 1 
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The following assumption and lemma are instrumental 
for our future developments. 

Assumption 4.2 The *-th row of factor T{q) is inde¬ 
pendent of q for i £ kPH 

Lemma 4.2 Under Assumption 4.2, there exists con¬ 
stant matrices Yj £ R nx,s , j £ n , such that, for all vec¬ 
tors z = col( 2 j) £ R" we have 

n 

Ru(q)z = {J2 Y jZj)r u - (11) 

j =i 

Proof 4.1 From (10) it follows that 

n 

D\ u = J2 e ieJ{ejGr u ). (12) 

2=1 

Using the definition of R u (q) we get 


R u (q)z = T T (q) 


Y,eieJ(eJCr u ) 

. 2=1 


T(q)z, 


= E TT (<l) e i e i T ( < l) Ze i Cr - 
2 = 1 


=E 


Y J T T {q)e i ejT(q) 

1=1 


ezj e,: Cr u 


(13) 


4-3 First robust momenta observer 

Proposition 1 Consider the system (1), (2) where the 
inertia matrix M{q) and the friction matrix 91 verify 
Assumptions 4.1 and 4.2. The 2n + s dimensional I&I 
adaptive momenta observer 

p I = -T T {q)[\/V-G{q)u-d} 

n 

- (E Y iPi)ru - [A Q(q) + Rk(q)]p 
2 = 1 
n 

fu, = Pi) {pi + Xp) 

2=1 

di = T{q)p 
P=Pi + A Q(q) 



2=1 


d = di + q 
P = T~ T (q)p 

with the constant n x n matrices Li given by (15), Q{q ) 
given in (9), Y, L £ R” xs , i £ n given in (14) and A > 0 a 
free parameter, ensures boundedness of all signals and 

lim [p(t) - p{t)} = 0. (16) 


for all initial conditions (q( 0), p(0)) £ 1" X R n . 


Hence, (11) follows swapping the sums and defining 


Proof 4.2 Let the observation and parameter estima¬ 
tion errors be defined as 


n 

Y j :=^2T T {q)e ie jT{q)e je JC 

2=1 

= E Lie 3 e JC, j £n (14) 

with matrices L,; defined as 

Li := T t { q)aejT{q) , i£n (15) 


It only remains to prove that the matrices Yj and Li are 
constant. Towards this end we refer to (14) and notice 
that, in view of Assumption 4.2, the term ejT(q) is con¬ 
stant for i £ k while the term ej C is an 1 x s zero vector 
for i k. Completing the proof. 

Remark 4.2 As indicated in Remark 2.1, except for 
the case when M (and, consequently, the factor T) are 
constant, to satisfy Assumption 4.2 we have to assume 
that some of the elements of r are known. See Section 6 
for some physical examples. 


3 Consequently, the unknown friction coefficients are located 
in these rows. 


p = p-p 

f u =r u - r u (17) 

d = d — d. (18) 

Following the I&I adaptive observer procedure [1] we 
propose to generate the estimates as the sum of a pro¬ 
portional and an integral term, that is, 

p = Pi + pp{q) 

f u = r UI + r up (p) (19) 

d = di+dp(q ), (20) 

where the mappings 

p P : R n -> R" 
r up : R” —>• R s 
d P : R" R n , 

and the observer states pi,di £ R n and r UI £ R s will be 
defined below Q 

4 The reason for the particular selection of the arguments 
of the proportional terms will become clear below. 


4 









First, we study the dynamic behavior of p and compute 

P = Pi + Vp P T(q)p + T t (g)[W - G(q)u] 

+ [Rk (q) + Ru(q)\(p - p) - T T (q)(d - d), 

where we have invoked Assumption 4.1 that ensures— 
via (6) and (8)—that J(q,p) = 0, and used (17) to obtain 
the terms in the second row. Invoking (11) we can write 


yields 

d = - Vd P T{q)p . (24) 

We will now analyze the stability of the error model 
(22), (23) and (24) with the aid of the proper Lyapunov 
function candidate 

V(p, d, f u ) = i(|p| 2 + \d\ 2 + \f u \ 2 ). (25) 


n 

Ru{q)p = (^2 Y iPi) r n- 

i =1 


Hence, proposing 

PI = -Vp P T(q)p - T T (q)[VV - G(q)u] 

n 

- f E YiPi) f u - Rk ( q)p + r T (q)d, (21) 

i= 1 


Taking its time-derivative we obtain 


V = -p 1 


R(q) + A i n P- P T r(r YiPiYu 


2=1 


- [ArJVr up +d r S7d P T( y q)]p. 


T T {q)d\ 
(26) 


Clearly, if the mappings r up (p) and dp(q) solve the par¬ 
tial differential equations (PDEs) 


yields 


P = ~[R(q) + V q p P T(q)]p- (£Y^)f u+ T T (q)d 

n 

= ~[R{q) + A In]p ~ (Y, Y iPi)?u + T T (q)d, (22) 
2=1 


= -jCYY Pi) 

2=1 

Vdp = I n . 


(27) 


one gets 


V = —p T [R(q) + XI n ]p < —A|p| 2 . (28) 


where to obtain the second equations we have selected 

Pp(q) = A Q{q), 

with Q(q) given in (9). 

Now, the time derivative of f u is given as 

f u = f UI + \7r up p 

= f UI + Vr„ p \pi + \7p P T(q)p\ 

= f UI + Vr„ p \pi + A (p - p)]. 


From (25), (28) we conclude that p G £ 2 fl £oo an d 
d , f u £ Coo ■ Doing some standard signal chasing it is 
straightforward to prove from here that (16) holds. 

Motivated by the conclusion above let us now study the 
PDEs (27). The second one has the trivial solution dp = 
q. Regarding the first one, it is clear that the elements 
of the mapping r up (p) must be of the quadratic form 

(r UP (p))i = ^fP TL iP’ i e 


Hence, choosing 


f UI = -Vr up (p/ + A p), 


yields 


f u = -A Vr up p. 


Finally, the time derivative of d is given as 


with constant, symmetric matrices Lj £ R nxra . Replac¬ 
ing the expression above in the PDE (27) yields 


(23) 


p T Li 


P T L S 


— 


d = di + X7dpT(q)p 
= d,! + S7d P T{q)(p-p). 


Hence, choosing 


That lead us to the solution 


-£ / i &j • • • L s Cj 


= ~ Y i > J e n. 


(29) 


di = ~'Vd P T(q)p, 


It only remains to show that the resulting matrices Li 


5 



are symmetric. From (29) we get 


whose solution is not obvious. 


Lj Lj ei Lj €2 • ■ • Lj e n 

Y\ ej 1 2 - * * Y n ej 

Clearly, the matrix Lj is symmetric if and only if 
ejYk = e^Yi, V?’, k G n, i ^ k. 

This fact can be easily verified using (14) 

n 

ejYi = el ^T T (q)ejejT(q)eieJC 
i=i 

n 

= e JYl TT (9)ejejT(g)e fc eJ C 
i =1 
= ejY k . 


5 ROBUST OBSERVER FOR GENERAL 
PERTURBED MECHANICAL SYSTEMS 
WITH KNOWN FRICTION 

In this section we redesign the I&I speed observer of [13], 
see also [2], to ensure its global convergence in spite of 
the presence of the unknown disturbances d and known 
friction forces in all coordinates. 

Proposition 2 Consider the system (1), (2) with 
known friction matrix 91. There exist smooth mappings 



such that the interconnection of (1), (2) with 


Replacing all the derivations above in pj, di and r UI 
gives the equations given in the proposition completing 
the proof. 

Remark 4.3 If Assumption 4.1 is not imposed a term 
J(q 1 p)p appears in the error equation (22) and (23). 
Even though this term is quadratic in the unknown state 
p, the properties of J{q,p) can be used to handle this 
term in the first error equations—this is done in the 
second observer in the next section. However, there is no 
obvious way to create a suitable error term for the second 
error equation styming the relaxation of Assumption 4.1. 


X A(X, q 1 u ) /01 \ 

p = B(X, q), 

where X € R 4 ™ x R>o, p G R”, ensures (16) holds for 
all initial conditions 

( 9 ( 0 ),p(0),X(0)) G 1" x IT x R 4 " x R> 0 . 

This implies that, in spite of the presence of the unknown 
disturbances d , (31) is a globally convergent momenta 
observer for the mechanical system with friction (1), (2). 


Remark 4.4 If the matrices Yj are not constant the 
first PDE in (27) does not admit a solution, hence As¬ 
sumption 4.2 is required. It can be shown that making 
r up function of q does not solve the problem, because a 
quadratic function of p will appear in V. 


Pro of 5.1 The construction of the observer follows very 
closely the one reported in [13] with the only difference 
of the inclusion of an adaptation law for the unknown 
disturbance parameters d. However, for the sake of com¬ 
pleteness, a detailed derivation of all the steps is given. 


Remark 4.5 As shown in Proposition 6 of [16] the dy¬ 
namics of mechanical systems satisfying Assumption 4.1 
expressed in the coordinates ( Q,p ) take the form 

Q = P 

P = -R(Q)p - f t (Q)[VV(Q) - G(Q)u + d], (30) 

where (• )(Q) := (Q)), with Q 1 : R" — > R” a left 

inverse of Q(q), that is, = z for all z G R”. 

Although the construction of an observer for (30) when 
the friction is known is straightforward, the case of un¬ 
known friction is far from trivial. Applying the I&I pro¬ 
cedure used in Proposition 1 leads, for the definition of 
dp(q), to a PDE of the form 

VS(Q) = T(Q), 


Define the estimation errors 

p = p-p 

d = d — d. 


(32) 


Following the I&I adaptive observer procedure [1] we 
propose to generate the estimates as 


p:= p! +p P (q,i, {>) 
d:= d/ + dp(q, r ) 


(33) 


where the mappings pp : R" x R" x R" —> R n and 
dp G R”, and the observer states d/ G R n , pi G R” and 
r G R are defined such that (16) holds. 
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We, therefore, study the dynamic behavior ofp and com¬ 
pute 


Now, since the term in brackets in (38) is equal to zero 
if £> = p and <4 = 9, there exist mappings 


P = Pi + ^qPpq + V 4 pp4 + 

- J(q,p)p + T r (q )[W - G(q)u] + R{q)p - T t { q)d. 

In [2] it has been shown that the mapping J{q,p) defined 
in (6) verifies the following properties: 

(P.i) J(q,p) is linear in the second argument, that is 

J{q,aip + a 2 p) = aiJ(g,p) + a 2 J(q,p) 


A q , A p : R" x 1" x -x K" xn 

verifying 

A,(«,t>,0) = 0, A p (q, 4>, 0) = 0, (39) 

and such that 

H(q,p) -H( 4, i>) = A q (q,i,e q ) + A p (q, 4>,e p ), (40) 


for all q, p, p € R ra , and aq, a 2 £ 1 
(P.ii) There exists a mapping J : 1" x R" 
ing 

J(q,p)p = J{q,p)p. 


satisfy- 


where 

e q :=i-q, e p := jp - p. (41) 

Substituting (36), (38) and (40) in (35), yields 


Hence, proposing 

Pi ■= -V 4 pp 4 - Vj,ppi) + J(q,p)P~ R{q)p- 
- T(q) T (q)\7V + v - \7 q p P T(q)p + T r (q)d, 

(34) 

together with Properties (P.i) and (P.ii) yields 

P = [J{q,p) + J(q,p) - R{q) ~ V g PpT{q)}p + T T {q)d. 

(35) 

It is clear that if the mapping pp solves the PDE 
Vqp P = [ipl n + J(q,p)}T~ 1 {q), 
with ip > 0 a design constant, the p-dynamics reduces to 

p = [J{q,p] - ipin - R{q)}p + T T {q)d. 

Recalling that J(q,p) is skew-symmetric and R(q) > 0 
the unperturbed part of the error dynamics above, i.e. 
when d = 0, is exponentially stable. 

Similarly to [13], to avoid the solution of the PDE, the 
dynamic scaling technique is used. Towards this end, 
define the mapping 

U{q,p) := [ipI n + •/(<?, p)]T _1 (g). (36) 

and define pp as 

pp(qA, i 3 ) ; = M-(A,i>)q- (37) 

The choice above yields V g pp = 'H( <4, 4>), which may be 
written as 

\7qPP = U(q,p) - [H(q,p) - U(A, i>)]. (38) 


p = lJ(q,p ) - ip in - R]P 

+ ^qiqA^q) + Ap(g, 4 pe p )^T{q)p + T T (q)d. 

The mappings A q , A p play the role of disturbances that 
are dominated with a dynamic scaling and a proper 
choice of the observer dynamics. For, define the dynam¬ 
ically scaled off -the-manifold coordinate 

V = -P, (42) 

r 

where r is a scaling factor to be defined. The dynamic 
behavior of q is given by 

f] = (J - R — I)rj + (Ag + A p )Trj+ -T T d- - 77 , (43) 

r r 

where, for brevity, the arguments of the mappings are 
omitted. 

Mimicking [2] select the dynamics of <4, 4> as 
4 = T(q)p-ip 1 e q 

4 = —T T (q)W + v + J{q,p)p - Rp (44) 
- ip 2 e P + T T (q)d 

where ipi , ip 2 are some positive functions of the state 
defined later. Using (44), together with (41), we get 

e q =T( y q)rir-ip 1 e q 
e-p = V q p P T(q)qr - ip 2 e p . 

Moreover, select the dynamics of r as 
r = -^(r-l) + ^(||A p T|| 2 + ||A 9 T|| 2 ), r(0) > 1, (46) 
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which clearly solves the PDE 


with || • || the matrix induced 2-norm. At this point we 
make the important observation that the set 

{r £ R : r > 1} 

is invariant for the dynamics (46). Hence, r(t) > 1, Vt > 

0. 

On other hand, taking the time-derivative of d, we get 

d = dj + V q dpT(q)p + \/ r dpf 
= d/ + \7 q d P (q,r)T(q)(p - rr[) + \/ r dpr 

and choosing 

dj = —V q dpT{q)p — \7 r dpf, (47) 

the d-dynamics take the form 

d =—rV q dpT(q)ri (48) 


i - rX7 q d P (q,r) = 0. 

It only remains to study the term rr, which is given by 

rf = -^r 2 +r^\\A pT \\ 2 + \\A qT \\ 2 ). 

Now, (39) ensures the existence of mappings A q , A p : 
I" x I" x R" —> R raxn such that 

\\A q {q, txejll < II A q (q, t>, e 9 )|| |e,| 
l|A p (g, t>.e p )|| < ||A p (g, i>,e p )|| |e p |. 

Hence 

||A P T|| 2 + ||A 9 T|| 2 < ||T|| 2 |(||A p || 2 |e p | 2 + |AJ 2 |e 9 | 2 ). 


We now analyze the error system (42), (45), (46), (48)— 
with the coordinate f = (r — 1). For, define the proper 
Lyapunov function candidate. 


V(r),e q ,e p ,r,z a ) := ^ (jry| 2 + |e g | 2 + |e p | 2 + r 2 + 
Taking its time-derivative we obtain 


Id?) • 

(49) 


V < ~ - 1^ \V I 2 - (Vi - ^r 2 \\T(q)\\' 2 ^j \e q \ 2 

- ^2 - ^ 2 ||V 9 pp|| 2 ||T(<7)|| 2 ^ |e p | 2 + fr+ 

+ d 7 ^ - rV q dp \ T(q)r/ 


Clearly, if we set 


Setting 


1p3 = K 


rr 


^4 = 777——imi1|A,f+ * 


^5 = 


4(1 + i>z) 

rr 

4(1 +^ 3 ) 


ni 2 iiA P n 2 + K , 


one gets 


V<-«(H 2 + |e,| 2 + |e p | 2 +r 2 ] (52) 

for some positive constant k. 

The proof is completed invoking the arguments of [13], 
selecting the observer state as 


= 4(1 +^ 3 ), il>i = \r 2 \\T(q)W 2 + ( 50 ) 

and 

V>2 = ^ 2 ||V (? pp|| 2 ||T(g)|| 2 + V’s, 

where V ; 31 V ; 4 , V ; 5 are positive functions of the state de¬ 
fined below, one gets 

V < -^ 3 |??| 2 - ipA\e q \ 2 - tp 5 |e P | 2 + rr 
+ d 1 ^ - rV q d P (q, r) j T(q)rj. 


X := (4, f>,pi,di,r — 1), 

and defining A(X, q,u) from (34), (44), (46) and (47) 
and setting B(X, q) via (33). 

Remark 5.1 It is clear from the proof that the key step 
to reject the disturbances is to make the proportional 
term of the parameter estimator, dp a function of the 
dynamic scaling factor r, see (51). 


6 PHYSICAL EXAMPLES 


To eliminate the cross term appearing in the last term 
of the right hand side above we select 

dp(q,r) = \q, (51) 


In this section we present three physical mechanical sys¬ 
tems that satisfy the conditions of Proposition 1. Con¬ 
sequently, robust adaptive speed observation is possible 
for them. 



6.1 Constant inertia matrix 


and the definition of all constants may be found in [16]. 
The Cholesky factorization is given as 


In the case of constant inertia matrix Assumption 4.1 
is trivially satisfied, because the factor T can be taken 
to be constant. Assumption 4.2 is also satisfied with 
dim(g) = n and C = I n , hence all friction coefficients 
can be identified. 

Given any constant factor T, the vector field Q(q) that 
solves (9) is given by 


Q{q) = T 1 q. 


Finally, from (14) and (15) we get 


Yj = ^2T T (q)eiejT(q)ejeJ j £ n 

i =1 


U = [CZt] T (T)\ * G n. 


6.2 Planar redundant manipulator with one elastic de¬ 
gree of freedom 

This is a 4-dof underactuated mechanical system de¬ 
picted in 1. The inverse inertia matrix is given by 



Fig. 1. Planar redundant manipulator with one elastic degree 
of freedom. 


M-frq) = 


7 ~T 0 0 

* ^ + 7 YH Cl2 

* * At 0 

a l 

1 


where we defined 


S 12 ■■= sin(<?i + q 2 ), C 12 := cos(gi + q 2 ) 

y/ Mml 


a 2 ■■= 


y/m + M 


, 03 := y/M + m, 


T(q) = 


1 

Vi 

0 

0 

0 

-1 

C~I 

J_ 

_ 02_ 

0 

0 

0 

P^--s 12 

/ m as 

J_ 

d3 

0 


-C 12 0 — 

m as as 


and the vector field Q(q) that solves (9) is 


Q(q) = 


Vlqi 

a 2 {qi + q 2 ) 

— 02 Y -^f<-U2 + a 3<?3 
/Ml c . 

-02V -±^12 + a 3 g 4 


A matrix C that satisfies Assumption 4.2 is 


C = 


I 2 

02x2 


Hence, we can consider as unknown the frictions in the 
elastic coordinate r± and the revolute joint r 2 . 


Finally, 


IV = 


Y 1 — 


u = 


\ 0 0 0 

T-^00 


0 0 0 0 
-fj «i 0 0 


7 0 


0 2 


0 0 

02x2 02x2 


Lo — 


1 
I 
a 2 
“VT 


CL2 

V7 

al 


02x2 

02x2 02x2 

6.3 2D-Spider crane gantry cart 


This is a 3-dof underactuated mechanical system de¬ 
picted in Fig. 2. The inertia matrix is 


M{q) = 


m r + m 0 mT 3 C 3 

* m r + m 171L3S3 

* * ?ttA§ 
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From the definition of T(q ) in (53) it is clear that 



Fig. 2. 2D-Spider crane gantry cart 


with inverse 


M~\q) = 


m r -\-mC^ 

(m r +m)m r 


171C3S3 

(m r -\-m)m r 

mr+m-mCg 

(m r -\-m)m r 


-c 3 

L/3m r 

_ 

m r L3 

m r -\-m 

m r Lim 


c T = [0 01 

satisfies Assumption 4.2. Hence, we can consider as un¬ 
known the friction parameter r 3 . 

Finally, from (14) and (15) we get 

n T = [0 0 c 2 ] , 

0 0 0 

Li = 0 0 0 

0 0 c 2 

We show several simulations of the proposed observer. 
The system has the two forces shown in Fig. 2 as con¬ 
trol inputs, that is, u = col (F x , F y ) and constant input 
matrix G of the form 


where, to simplify the notation, we have defined 

S 3 := sin(g 3 ),C 3 := cos(g 3 ), 

and the definition of all constants may be found in [9]. 
An upper triangular factorization of M~ 1 {q) is given as 

a 0 — bC 3 

T(q) = 0 a -bS 3 , (53) 

0 0 c 

where we defined the constants 



We can check that the columns of T(q) satisfy ( 8 ) and 
thus the system verifies Assumption 1. Taking the in¬ 
verse of T(q) we get 

\ 0 aL 3 mC 3 

T~ 1 {q) = 0 i aL 3 mS 3 , 

OO i 

From the equation above it is clear that a mapping Q(q) 
that solves (9) is 

\qi+ aL 3 mS 3 
Q(q) = \q -2 - aL 3 mC 3 ■ 

c«3 



The parameters are taken as m r = 0.5 kg for ring mass, 
to = 1 kg for payload mass and L 3 = 0.5 to for the cable 
length. We fix the control inputs as F x = 1.535 cos(t) 
and F y = 7.67sin(t). The disturbances are taken as 
d = col( 0 . 1 , 0 . 2 , 0 . 2 ) and the friction coefficients r = 
col(0, 0, 0.5), with r 3 being unknown. 

The transient behavior of the error signals p, f 3 and d 
with the tuning parameter A = 0.8 and different initial 
conditions of di , r 3l and Pi are shown in Fig.3 and Fig.4. 
As seen from the figures, besides the convergence to zero 
of the momenta estimate predicted by the theory, we also 
observe that the estimated parameters converge to their 
true value—assessing the fact that the signals chosen for 
the simulation are persistently exciting. 

To evaluate the effect of the tuning parameter A on the 
transient behavior we also show in Fig.5 and Fig .6 the 
transient behavior of the error signals for different values 
of A. 

Finally to illustrate the robustness of the adaptive ob¬ 
server, we carried out a simulation considering that the 
input disturbance d\ is subject to step changes. The tra¬ 
jectories of di(t) and its estimation, depicted in Fig. 7, 
clearly illustrate the tracking capability of the proposed 
observer. 

Remark 6.1 It is interesting to note that the standard 
(lower triangular) Cholesky factorization of M~ l {q) 
does not satisfy ( 8 ). 
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Fig. 6. Transient behavior of d and for different values of A 


Fig. 3. Transient behavior of p for A = 0.8 and different 
initial conditions of di , and Pi 



Fig. 4. Transient behavior of d and r 3 for A = 0.8 and differ¬ 
ent initial conditions of di, r^j and Pi 



Fig. 5. Transient behavior of p for different values of A 

7 Future Research 

The design of the observer in Proposition 1 requires the 
explicit solution of the PDE (9). This requirement re- 



Fig. 7. Transient behavior of d\ and di with A = 2 

stricts the practical applicability of the approach. In¬ 
deed, this PDE has no free parameters and its explicit 
solution may be even impossible. As indicated in [16] 
this is the case of the classical cart-pendulum example. 

Current research is under way to extend the realm of ap¬ 
plication of the observer in Proposition 1. In particular, 
it is possible to consider the following generalization of 
the Lyapunov function candidate (25) 


lif 


fJp~ 


\d\ 2 , 


with P > 0 a constant matrix. It is straightforward to 
show that the second PDE in (27)—that cancels the 
cross term appearing in the derivative of the new Lya¬ 
punov function—becomes 


Vr„ P =-t(E PF ^)’ 


A 


where we underscore the presence of the matrix P in 
front of Yi. There are inertia matrices where the corre¬ 
sponding Yi are not constant but there exists positive 
definite P that will make PY .] T constant- hence relaxing 
Assumption 4.2. We are currently investigating whether 
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there exist physical systems for which such property 
holds. 

Another, quite challenging, task is the extension of 
Proposition 1 to systems that do not have ZRS. One 
possibility is to look into the next class of systems par¬ 
tially linearizable via coordinate changes characterized 
in [16]. 
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